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We consider various one-dimensional non-equilibrium models, namely the diffusion-limited pair- 
annihilation and creation model (DPAC) and its unbiased version (the Lushnikov's model), the 
DPAC model with particle injection (DPACI), as well as (biased) diffusion- limited coagulation model 
(DC). We study the DPAC model using an approach based on a duality transformation and the 
generating function of the dual model. We are able to compute exactly the density and correlation 
functions in the general case with arbitrary initial states. Further, we assume that a source injects 
particles in the system. Solving, via the duality transformation, the equations of motions of the 
density and the non-instantaneous two-point correlation functions, we see how the source affects 
the dynamics. Finally we extend the previous results to the DC model with help of a similarity 
transformation. 
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I. INTRODUCTION 



Stochastic reaction-diffusion models play an important role in the description of classical interacting many-particle 
nonequilibrium systems in physics and in interdisciplinary areas (chemistry, biology, economics, ...). Usually physical 
systems are much too complex to be treatable analytically or even numerically. However, in the context of critical 
phenomena simple toy models have been shown to be useful in determining universal properties and understanding 
possible relationships between microscopically different processes (see e.g. |l[ and references therein). 

In this work we study, on a periodic chain, models which are prototypes of one-dimensional diffusion-limited 
reactions. Our main purpose is to present a novel approach (based on a duality transformation combined with 
generating function techniques [^|) to analyze the DPAC and related models. In this paper we illustrate this approach 
by recovering known results and deriving, in a simple and systematic manner, some new ones and show that the 
method employed here is complementary to the traditional (see e.g. [^f]?! and references therein) ones. We also solve 
the dynamics of the dual of the DPAC model which is a biased generalization of Glauber's model. We believe that 
the method used here is particularly suited to site-depend and/or disordered systems. A more detailed presentation 
of the method, as well as other applications, will be given elsewhere 

Diffusion-annihilation and coagulation models (in their free-fermionic version) properly describe the kinet- 
ics of excitons in several materials: the dynamics of photoexcited solitons and polarons in, e.g., chains of 
[Pt{en)2\[Pt{en)2Cl2]{BF4)i, called (MX), where (en) denotes ethylenedy amine, or the fusion of photogcnerated 
excitons in chains of tetramethylammonium manganese trichloride (TMMC) Q. Models of pair-annihilation and 
creation are useful to decribe problems of dimer adsorption and desorption H . 

The motivation for studying such one-dimensional systems is not only their experimental relevance, but also their 
theoretical importance in the understanding of the fluctuations of low-dimensional systems. In both one- and two- 
dimensional systems, diffusive mixing is inefRcient and leads to the building of large-scale correlations. The mean-field 
is not adequate and in this sense exact results are desirable. 

The paper is organized as follows: In the sequel of this section we recall the stochastic Hamiltonian formulation 
of Markovian processes obeying a master equation. In section II we introduce two models: the diffusion-limited 
pair- annihilation- creation (DPAC) model and the diffusion-limited pair- annihilation- creation model connected with 
a source (DPACI model). In section III we map through a duality transformation the DPAC and DPACI models on 
other stochastic models. The section IV and V are devoted to the detailed study of the DPAC model. In section IV, 
we evaluate the exact generating function of the dual DPAC model. Correlation functions of the DPAC model are 
studied in section V. In section VI, we study the DPACI model via the duality transformation. In section VII we show 
how to extend the solution obtained for the diffusion-limited pair- annihilation (DPA) model to the diffusion-limited 
(DC) model. The section VIII is devoted to the conclusion. 
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In this work we study some one-dimensional two-states nonequilibrium systems for which the dynamics takes place 
on a periodic chain with L (even) sites. The dynamics of particles (of a single species) is governed by a Master 
equation. Each site of the lattice can be either empty or occupied by a particle at most (the hard-core constraint), 
say, of species A . It is known that such nonequilibrium problems can be reformulated as field-theoretical many-body 
problems. Below we briefly recall the basis of the field-theoretical formulation of the master equation. 

The state of the system is represented by the ket \P{t)) ~ X]{ri} -^({"^i' where the sum runs over the 2^ 
configurations. At site i the local state is specified by the ket \ni) = (1 0)-'" if the site i is empty, \ni) — (0 1)"^ if 
the site i is occupied by a particle of type A. We define the left vacuum (x|, where {x\ = The master 

equation can be rewritten formally as an imaginary-time Schrodinger equation: ^|P(i)) — —H\P{t)), where H is the 
Stochastic Hamiltonian which governs the dynamics of the system . In general, it is neither hermitian nor normal. Its 
construction from the master equation is a well established procedure (see e.g. [D and references therein). Here we 
specifically focus on two-state systems so that the stochastic Hamiltonian can be rewritten in term of Pauli's operator 
as a (pseudo-)spin chain (see below). Probability conservation (stochasticity of H)yields{x\H — 0. 



II. THE DPAC AND DPACI MODELS 



In this section we present two models which we will consider in the sequel of the work. We begin with the diffusion- 
limited pair- annihilation- creation (DPAC) model. 

We consider, without loss of generality, a periodic chain of L (even) sites on which an even number of particles with 
hard-core repulsion move according to a Master equation Q 

In the DPAC model, at each time-step four events can occur: 

A particle can jump from site r to r -I- 1 (provided the latter was previously empty) with a probability rate h' . 
A particle can jump from site r to r — 1 (provided the latter was previously empty) with a probability rate h. 
A pair of particles can be created at sites r and r -\-\ (provided the latter were previously empty) with a probability 
rate e'. 

A pair of particles can be annihilated at sites r and r + 1 (provided the latter were previously occupied) with a 
probability rate e. 

The Master equation associated to these processes can be formulated as an imaginary time Schrodinger equation. 
Identifying a vacancy with a (pseudo-) spin up and a particle with a (pseudo-) spin down the stochastic Hamiltonian 
can be written in term of Pauli (pseudo-)spin— 1/2 operators, as follows: 

L 

ttDPAC \ " ttDPAC 

r=l 



rrDPAC I e'-e+{h-h') \ / e'-e-{h-h') \ f e' + e + h + h' . 



e' + e - h - h'\ . ^ ^ „ „ , . f e' — e — (h ~ h')\ „ „ . f e' ~ e + (h ~ h')\ „ „ 



where aif , trf are the usual Pauli's matrices acting on site j and 2af = a^±ia^ . The operator n, = i (l — erf) is the 
local (at site j) density operator (1 denotes the identity operator). This model (with a particle-hole transformation) 



was introduced by Grynberg et al. and corresponds to a biased generalization of Lushnikov's model 10 
The constraint 

e + e' = /i + /i', (2) 



'^Because the DPAC model preserves the parity of the number of particles, we should split the dynamics in a even and an odd 
sector. Here we focus specifically, and without loss of generality, on the even sector. The treatment of the odd sector is similar 
excepted for the boundary conditions (see e.g. [^). The parity of the number of particles does not play an important role for 
the dynamics of the model. The restriction to the even sector is performed projecting out the odd sector from the initial state, 
say |P(0)), with help of the projector | (1 + Q), with Q — cr^. In the sequel, for simplicity, we will often call |P(0)) its 
even-sector projection, i.e. |P(0)>'''"=" = \{1 + Q)|P(0)). 
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is the free-femionic condition, for which the model becomes soluble q In the sequel, we assume that the condition (Q) is 
fulfilled. When, in ([^), no pairs are created, i.e. e' = 0, the model under consideration is called DP A {diffusion-limited 
pair- annihilation) model. 

Let us now consider the following stochastic Hamiltonian (with periodic boundary conditions): 

L L 
r— 1 r—1 

This stochatic Hamiltonian term corresponds to a single particles "source" which injects in the system, with rate 
7, particles at site j < L and L, whenever both sites were previously vacant. If both sites j < L and L were 
previously occupied, the "source" annihilates (with rate 7) both particles. When one of the sites j < L or L is 
occupied and the other empty, the effect of (with rate 7) is to fill in the previously vacant site and to empty 

the previously occupied one. In (^), the term in (t£ is an artefact due to the periodicity of the problem (because 
of the duality transformation, see next section) and preserves the parity of the system described by the stochastic 
Hamiltonian ijDPACi = jjDPAC _|_ jjmput rpj^^ steady-state as well as the density relaxation time of such a system 
(in its translationally invariant version) have been studied for the case where there was input of particles (i.e. 
7 > 0), but no pair-creation (i.e. e' = 0), for an infinite system (where there is no problem of boundary conditions and 
one can consider simply _iJ™P"* = 7^^ [a^ — 1]). Here we obtain the exact density of the DPACI model and study 
the interplay between the pair-creation and the source (see (|3^)). We also obtain the two-point non-instantaneous 
correlation functions of the (unbiased) DPACI model in the absence as well in the presence of pair-creation. 



III. THE DUALITY TRANSFORMATION 

Following 1^,^, we introduce a set of operators which forms a periodic Temperley-Lieb algebra (a quotient of a 
Hecke algebra). Further, we define an unitary transformation which allows to map the DPAC model on another 
stochastic model for which we are able to calculate the generating function. 

Let us define the following operators in term of Pauli spin-matrices: 

62,-1 ^ ^(1 + ^J), 1 < J < i ; = ^(1 + a|a;+i), 1 < J < L - 1 ; eo = = ^(1 + Calat), (4) 

where (7 = nj=i...L<^f 

We define the (unitary) duality transformation V by 




2L-1 



F^exp -^aj l[[{l + ^)e,-l]\ (5) 



The dual (unitary) transformation V is also called domain-wall transformation [Qjj] because of its physical inter- 
pretation. In this duality transformation has been introduced to study the zero-temperature Ising model with help 
of its free-fermionic dual counterpart (the DPA model, where e' = 0). The transformation (||) connects a stochastic 
model with periodic boundary conditions to another stochastic models with the same boundary conditions (2) . 

Let us define the dual version of the stochastic Hamiltonian (|^): 

^DPAC ^ v-'H^f^f{e, e', h')V ; i?^^-^*^ ^ ^r^^^ (6) 
with periodic boundary conditions: crf^^j. = (jf, {r < L) 



^Recently we have proposed jl2| an approach to study analytically the DPAC model beyond the free-fermion case (g) . Although 
the approach provides the correct long-time behaviour of the density and correlation fuctions (including the subdominant 
terms), it turns out to be less general than expected |l3j. It is not the purpose of this work to study the limits of the approach 
devised in [Q. 
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From now on, we work with the dual model (H) and obtain information on the original model (with the constraint 
@)). Notice that the dual model (||) is still a stochastic one because, by construction, = {x\HJ^^'^'~^ = 

{x\H^^^'^' = {x\Hr ~ 0, since {x\V = (x| = j where a — — e'T^-^^-^^ @ is a constant which plays no 

role in the following and therefore will be omitted. 

Here we consider the domain-wall duality transformation and, formally, obtain (when (|^) is fulfilled) a generating 
functions for the dual model. The latter will allow to solve the DPAC model completely. In the third part of this 
work, we translate our results into the language of the diffusion-limited coagulation model. 

The duality transformation ^ maps the DPAC model (^ onto a biased generalization of the Glauber- Ising model 
1^. In the absence of bias, model (lb reduces to Lushnikov's model and the the duality transformation (^ maps this 
model onto Glauber-Ising model p]7|, p^ , p^ . 

Let us also consider the dual transformation of i/™P"* (||): 

L L 

jjtnput ^ J2 = 7^ (1 - . . . , (7) 

r=l r— 1 

with periodic boundary conditions, and thus H^paci ^ JjDPAC _^ fjinput _ 



IV. THE GENERATING FUNCTION OF (THE DUAL OF) THE DPAC MODEL 

In the first part of the paper we study exact properties of DPAC model using generating function of the dual model 
(^, which are explicitly computed in this section. We recover known results and produce new ones. 
We introduce Grassmann numbers "qm with their usual anticommuting properties: {?7m,?7ri} — 0, ym,n 
Following ||], we consider the quantity: 

G±({77, j},<) ^ (^1 n (1 ± V,<^]) e-^°"'"*|P(0)) = ( n (1 ± V,<^j) W ) ' (8) 
j>i \j>i I 

where \P(fi)) symbolically denotes the initial state in the original (Q) model with an even number of particles (see 
footnote 1), and |P(0)) = V~^\P{Q)) . G^{{r],j},t) is the generating function of the dual model and its derivatives 
provide correlation function of the dual model (||), e.g {afaj){t) — ^ ^gj^^gjj''.^'*^ I{'?}=0; (* < j)- 

It has to be emphasized that at this stage all the correlation functions obtained from the generating function 
G^{{ri,j},t) are correlators of the dual model (^). For the sequel, we introduce the following notations: 

h + h' + e + e' h - h' + e - e' , h' - h + e - e' ^ , , ~ , , , , 

= ^ ; c = ^ ; a = ^ \c— c + a~e — e]d=c— d = ri — ri (9) 

It is useful to separate the generating function into two parts. The one that generates the correlators with an even 
numbers operators cr^ denoted V"'"({?7, j}, i). The functional, called V~({?7, j}, i), generates correlators with an odd 
number of operators ct^, i.e., 

V± ({77, j}, t)^\ (G+ ({r;, j}, t) ± G- ({,7, j}, t)) (10) 

When the free-fermionic constraint (||) is fulfilled, the equation of motion of the generating function can be 
rewritten, using the properties of Grassmann algebra, as a first-order partial differential equation which can be solved 
by the method of characteristics : 



V^i{v,j},t) - ^ ( n (1 + - 0)) ± n (1 - = 0)) I ) exp 

\ \j>l j>l 



. ki ki>k2 



(11) 



where 77^ = e X]j^i(^^(^))/c/^?j (^^ thermodynamics hmit, L 00), 

(^*)jTfe(0 = ^y d0exp(i(j-fc)0+|2ccos0-de'*|ij =,5^Vfc(Ai) (12) 



4 



where the I^('w)'s are the usual modified Bessel function of first kind, z = — , A = \/cP ~ (P and S = f ^S^-^ ■ 

According to the definition of the generating functions (|8|),([fOD the instantaneous correlation functions of the dual 
model (^ are obtained by taking partial derivatives of V* {{rj, j}, t) (pi]). If one considers the mean- value of an observ- 
able O of the original DPAC process (§), its dual counterpart is denoted by {d{t)) as: {d{t)) = (x|6e"-^"^'*^*|P(0)) = 
{x\y{V'^OV){V-'^e-"'"'^''*V)V-^^\P{Q)) = (0{t)), and therefore, the dual counterpart of the observable O (say 



n 



, = i(l - a|)) is O ^ V-'OV (n, ^ i(l - a^_^&^)). 
For the unbiased DPAC model (i.e. Lushnikov's model, where h = h'), we thus recover Glauber's original result p6[ . 
For the biased DPAC situation, the generating function (|ll|) not only provides the complete solution of the DPAC 
model but also the complete solution of the biased generalization of Glauber's model with transition probability 
uj{ai —<Ti) = D {1 — -^a-iicri-i + ^Ci+i)), where the ai = ±1 are the usual spin variables and we assume the rate 

D = 6/2. For this model we obtain: (crj(t)) = ^^'j^J'^" l{.)}^o = e-''*J2ki'^k{0))S'^Ik-j{^t). Relabelling the sites 

which now run from + 1 to -j, the long-time behaviour follows as {(Jj{t)) — {a{0))e^^''~'^'>* ~ |--^^-^/((5), where 



d ^ = ^ = Vc^ ~ (P and f{u) = X^fe ^ ("'(O)))' which is considered to be an analytic function. 

This result generalizes recent results [Oj obtained for the biased zero-temperature Glauber's model. 



V. DENSITY AND CORRELATION FUNCTIONS OF THE DPAC MODEL 
A. Density and Correlation functions for arbitrary initial states 



In the previous section, we have obtained an explicit expression for the generating function of the dual model (see 
(^) and we have shown how to compute correlation functions for the dual model (^. In this section we show how to 
relate correlation functions of the dual model (||) to the correlation function of the orig inal DPAC model (§). 

Here we are especially interested in density-density correlation functions: {rij^ . . ■nj^){t) = ^((1 — tT|^_i(TU(l — 
'^h-i'^h^ ■ • • (1 — (ji < J2 < • ■ ■ < jn), where the symbol means that the mean-value {(T^){t) is 

taken with respect to the dual model (^. 

In what follows, we need to know how to connect initial correlation functions of the dual model (0) with initial 

those of the original DPAC model @): {a]^ . . . <t|^J(0) 
In particular, for the density, we have 



(n,i (1 - 2n,(0))) , {31 < J2 < . . . < J2„) 



K)(i) 



1 - 



1- z 



~2bt 



/„-i(2a)-/„+i(2a)} 



n>0 



-2bt 



E ( n (l-2n,(0))U 

J2>il \jl<j<j2 I 



jl+J2-2r+l 



{/,,+i_,(At)/j,_,(At) - /,,_,(At)/,,+i_,(Ai)} , 



(14) 



where A = \/(? ~ (P, S ^ = ( ^-i] and we have relabelled the sites indices which now run from — + 1 to . 

\c+d/ ... 22 

We can also obtain the non-instantaneous two-point correlation functions : 

1 



{nr{t)n,iO)) = 



1- z 



-2bt 



J2 { In-li2ct) - In+l{2ct)} 



n>0 



-2bt 



E ("^(0) n (1 - 2n,(0)) ) 5^^^^""^' {/,,+i_,(Ai)/,,_,.(At) - /,,_,(Ai)/,,+i_.(At)} (15) 



l2>ll 



jl<j<j2 



To obtain the long-time behaviour of (|l5|), we introduce the Fourier transform of the initial state from which we have 
substracted the homogeneous part: 

/(u, «, s) = e'("^'i+"J^) ((n,(0) Uj,<j<j, (1 - 2n,(0))) - (1 - 2^(0))^"^--'"^) Q{,n ~ Ji), where e(j2 - Ji) = 1 if 

j2 > ji and vanish otherwise. With the steepest- descendent method, we obtain for |r|, \s\ <^ L 00 and bt, |c|<. At 3> 1: 
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{nr{t)ns{0)) ~ p{t) 



-2{b-A)t 



(16) 



(S-'-(Ai)2- 

where p{t) is the translationally invariant and uncorrelated density studied in (17-25|). From ([l^), we see that the 
decay of the connected correlation functions (^5|) depends on the bias and is in the form of an exponential, except 
for the unbiased case (i.e. d — 0) and when there is no pair-creation (i.e. 6 = c). In the latter case the decay follows 
the power-law ~ (ci)^^. For the case d < 0, because of the drift, the rightmost sites tend more rapidly to their 
steady-state. The effect of the initial state appears through the function /, which is assumed to be analytic. 
All the the multipoint correlation functions can be obtained in a similar and systematic way. 



B. Density for translationally invariant and uncorrelated initial states 

In this subsection we check the results for the density (^4|) in the translationally invariant DPAC model against 
previous results obtained for some particular initial states (lattice with initial density p(0) = 0, 1/2, 1), directly by the 
free-fermionic approach P jTl| jl5| , p0| , ^ . As original results, we obtain the exact asymptotic behaviour of the density 
for uncorrelated and translationally invariant states of arbitrary initial density (in the presence as well as in the 
absence of pair-creation) . 

As one can check by direct computation, for tranlationally invariant systems, the density in the DPAC model is 
independent of the bias d and thus the density of Lushnikov's model coincides with the density of the generalized 
DPAC model (for the class of translationally invariant systems under consideration here). It has been shown to be a 
general property of the DPAC model, which holds true for every instantaneous correlation functions [T^ . We have 
seen in ( p^ , p^ ) that this property is lost for non-instantaneous correlation functions. 

For translationally invariant systems and with help of the properties of Bessel functions, (ITj) leads to {d'^<T^,^){t) = 



(tTgaJ)(<) ^z' + e-^"" Es'>o (('^0^.' )(0) - ^"') {Is-A^ct) - h+A^ct)), {s > 0). Where, z = ^-^Jf^\ This result 
coincides with Glauber's original one p^ . 

Let us consider the uncorrelated case where initially {ni^ . . . ni^){0) — (p(0))", (in > in-i > ■ ■ ■ > h)- This implies: 
= ((1 - 2rii) ... (1 - 2ns')){Q) = (1 - 2p(0))^' and therefore. 

Pit) - 2 (1 - {<^'o^l){t)) = - E ((1 - 2p(0))^ - ) (/..-i(2a) - /..+i(2a)) (17) 

s'>0 

To study expression ( p^ ) we notice that it can be rewritten as 

Pit) - p{^) = £1!!^ J [d^i _ y)i,{2m + iz' y')h{2m] + e ((i - y')^"'' - (i - ^')\^r') um) \ m 

[ n>2 J 

where y = (1 — 2p(0))signc. 

For the case where c — with e = e' > the density (|l^ ) reads p{t) — ^ + ip{0) ~ ^)e~^^*. From now we focus on 
the cases where |c| 7^ 0. 

In the absence of pair-creation (i.e. e' = 0, & = c = e > 0), we simply have (imposing z = 1, in (^8|)) 



Pit) = piO)e 



-2ct 



Ioi2ct) + - 2p(0))"-i/„(2a) 



(19) 



this result coincides with the result obtained in |19[ |. 

Taking into account the asymptotic behavior of the Bessel functions /„(x) and collecting terms, we arrive at the 
asymptotic behaviour of the density. 

In the massive (non-universal regime) when bt > \c\t ^ l,u — L"^ /\c\t ^ 1 and z < 1 < p(0) < 1, we have: 



P\t) ~ P(oo) ^ , 

32|Sl<y^ 



n I ^1 , |z|(|z| + l)(15-8|z|-3z2) v{v + l){lb-8y~iy^) 
i\z\ - y)il - 3(2/ + |z|)) + 1^1)2 (73^)2 



(20) 



The validity of ( pO| ) is restricted to \y\ < 1, i.e. to < p(0) < 1, which corresponds to the convergence domain of 
the geometrical series occuring in computing (pO|). The cases p(0) = and p(0) = 1 correspond to y = 1 (for p(0) = 0) 
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and y = —1 (for p(0) = 1). As \z\ < 1, the residual summation over the expansion of z" ^/„(2|c|i) can be carried 
out, leading to 



p{t) - p{co) 



I iv\c\t 



, if p(0) = and ? > 



32( 

-2(b-|<:|)t 



(l - 3(|z| - 1) + MiH^(i5 - 8|z| - 3z2)) , if p(0) = and 
1 - 3(|z| - 1) + M^U^(15 - 8|z| - 3^2)) , if p(0) = 1 and 



32(7r|c|t23/2 

„-2(t-|c|)f 



(21) 



/7r|c|t 



if p(0) = 1 and H < 



Results (^ ) coincide with results obtained in 

In the critical case, when there is no pair-creation (e' = 0, 5 = c > 0, z = 1), we have for bt ~ ct ^ l,u = L'^ /\c\t <C 1 
and \y\<l^ p{0) < 1. 



Pit) 



1 



2Vnct 



16ct 



(l-y)(3y + 2)-y(l + 2/) 



+ o ((a)-5/2) 



(22) 



This result is restricted to < p(0) < 1. When initially there are no particles, p(0) — Q,y = 1, no dynamics take 

(i + o((a)-i)). 



place. The case of an initially full lattice {y = — l,p(0) = 1) yields: /5(i)|p(o)^i = e *-^o(2c<) — 

A similar asymptotic results would have been obtained for the case where there are only pairs created and no 
annihilation (5 = c = e' > and e = 0). In this case z = —1 and p(oo) = 1. With (p8|), we obtain in the 
asymptotic regime (e't ^ 1, for an initially partially filled lattice, < p{Q) < 1) a critical decay of the density: 



Pit) 



1 



1 

16bt 



+ 0((6t)-5/2). 



_(l-y)(3y + 2)-y(l + y)M(«^ 

So far we have considered asymptotic behaviour for times which were much larger than the typical times of diffusion 
(we had u = 2L'^/\c\t ^ 1). Now, we study the asymptotic behaviour of the density, both in the massive and critical 
regimes, for typical times of order of the diffusion time, i.e. u — 2L'^/\c\t ~ 1. 

In the massive, non-universal regime, when bt > ct ^ l,u = 2L'^/\c\t ~ 1, from (p^, we obtain (Vy, i.e. VO < 
PiO) < 1): 



P{t) - p(oo) 



-2{b-\c\)t 



\z\i\z\ + l)-y{y + l) + Y,i {(1 - y2)y"-i _ (1 _ z^)\zr^}e~ + O {{Wr') 



In the critical regime when e' = and bt — \c\t 3> 1, with u = L"^ /\'c\t ^ 1, from (Eq), we obtain (0 < p(0) < 1): 



(23) 



Pit) 



1 



2V7rct 



1 - 



y(y + i) 



n>2 



O 



(a)-3/2 



(24) 



If one considers 



< 



= u/Stt <C 1, we recover the universal regime (p2[). A similar study can be performed 

e' > 0, for a partially filled 



47r|c|t 47r|c|t 

for the case where there is no pair- annihilation, but only pair-creation (e = 0, 6 = c 



initial lattice: p{0) < 1) and, we have p{t) ~ 1 — 



1 - 



kzvl 



, y" 



2\ TTCt 

Let us conclude this subsection with the study of the so-called staggered current for the (bias-)DPAC model. This 
quantity has been introduced by Grynberg et al. to measure the flux of particles due to the bias (so it vanishes for 
Lushnikov's model): it is defined as S{t) = {h'rimi^ — Tim+i) — hnm+i{i — nm))it). In 1^ this quantity was computed 
for an initially empty lattice. Here we obtain the exact expression of this quantity for translationally invariant (and 
uncorrelated, but random) initial states with arbitrary initial density /o(0). 

For a translational invariant system, the expression of the staggered current is Sf{t) = —d{p{t) — (nm+i^-m) (^)) — 
— ((t| 0-1^2)) • With help of results of this section, we then find: 



Jit) 



p(oo)(p(^)-l) + e-2^*^ [(l-2p(0))^' 



s'>0 



(/,'_2(2a) 



is'+2 



(2a)) 



(25) 



The long-time behavior of this quantity follows similarly as in (|20|-[24| 
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C. Density for translationally invariant but correlated initial states 



In the previous subsection, we have obtained exact results for the ordered DPAC model for a class of translationally 
invariant and uncorrelated initial state. It is also instructive to consider translationally invariant but correlated initial 
states of the form J^({crjcrj_,_^})|0), which according to (||) are transformed into V~^T{{ajaj_^_i})V\0) = ^({crj})|0), 
where J-{0) is a functional of the operator O. For translationally invariant systems the density is independent of the 
bias In the sequel, we focus on two classes of translationally invariant but correlated initial states of the form 
•^({^;^J+i})|0), namely: 



i) 



\Pm ^ n ( ^ + V^->J+i ) |0) (26) 



We have the initial correlations: ((7|)(0) — /i, (cr|cr|,_^^) (0) ~ and so the states ( p6| ) correspond to uncorrelated 
states for the dual model. However, via the duality transformation , this state is related to a correlated initial state 
in the original DPAC model (|]) QJ^. For the density of the DPAC model, with initial states (p6|), we have (see (IlT 



p{t) = + —-Y^z^ {Is-i{2Zt) - /,+i(2a)) - (/o(2ct) + h{2Zt)) , (27) 

s>0 

When c = 0, with e = e' > 0, this expression reduces to p{t) — ^ . Hereafter, we focus on the more 

interesting case c 7^ 0. 

In the first term of the r.h.s of (p7|), we recognize the expression of pit) — p{oo) for an initially uncorrelated state 
with p(0) = 1/2 (|l7|). When there is no creation of pairs of particles (i.e. e' = ^ 6 = c = e > 0, z = 1), 

P{t) = [I, {2ct) + h (2ct)] (28) 

In fact, such initial {correlated) states have been considered previously by Santos ||^ who computed the density for 
the diffusion-annihilation version of the DPAC-model (i.e. e' = ^ z = 1) using Jordan- Wigner transformation and 
the free-fermionic procedures. 

Let us study the asymptotic behaviour of the density for this (correlated) initial state (|2^). To this end, we 
proceed as in the previous subsection: we begin with the massive, non-universal regime, where bt > \c\t 3> 1, with 

u — /\c\t <^ 1. In this regime, the density decays as p{t) — p{oo) — — 7==- (m^(1 + signc) + 0{{\c\t)~^)) . On the 

4\/ TT\c\t 

other hand, in the critical regime (when e' = 0, 5 = c = e > 0,S 3> 1, with u ~ 2i^/|c|i <C 1), we have the following 
power-law decay: p{t) ~ -^—7^ ( 1 ) 

2^/ -net V Set/ 

As noticed in [g| , we see that although the initial state is correlated, the long-time behaviour of the density decays 
as t~2, as in the uncorrelated cases. The interesting point for the DPA model (where e' = 0) is however that the 
dynamics, though critical, is no longer universal: the amplitude of the density (i.e. the term 1 — p^) depends on the 
initial state. We infer that in this case initial correlations don't renormalize the dynamical exponent, yet affect the 
amplitude in a non- universal manner. 

Similarly, a power-law decay is obtained in the case where there is no pair-annihilation (i.e. b = —c = e' > 0): 
p{t) '--^ 1 — . Notice however that the asymptotic behavior is independent of p and so is universal. 

ii) Let us now consider another class of translationally- invariant, but correlated, initial states, namely (p > 0): 

for which initial correlation functions are: (CTf^nO'n)(0) — e~'^ sinh/3 ( 1 -I- %i -I- (1 — (5, n — 2(5, 1 )e^'^ sinh/3 ) and 



<^f>o^o>(0) = e-'5sinh/3(l + %i + (l-,5,,o-2<5,^i)e''3sinh/3) 



((t|)(0) = e ''sinh/3, where /? = In^. The initial states are therefore correlated both for the dual model and 
for the DPAC model. 

Notice that here the operators coding correlations act are not nearest-neighbours preventing a direct Jordan- Wigner 
transformation of the expression (|2^), i.e. a direct free-fermion approach. 
The computation of the density yields: 



8 



1 _ g-26t -26t-/3 

p{t) = ^— + ^— ^ {Is-i{2ct) - Is+i{2ct)) sinh/3(l + e-" smh/3) {Io{2ct) + /i(2a)) 



2 2 

s>0 



^ {2e-^ smh/3- l) e'^^* (/o(2a) - 72(2^)) (30) 



In the first term of the r.h.s. of (|3C|), we recognize the expression for the uncorrelated density p{t) ~ p{oo) with an 
uncorrelated initial state with p(0) = 1/2 (p7|). 

When c — with e = e' > 0, (jsy) reduces to the following expression: p{t) ~ ^ — e^'^sinh/3(| — e'' sinh/3)e~'*'^*. 
Below, we focus on the case where c ^ 0. 

When there are no pairs created (i.e. e' = 0, 6 = c = e>0 and thus z = 1), this expression reduces to: 

pit) ^[l-e^ sinh/3(l + e~'^ sinh/?)] — - {Ia{2ct) + h{2Zt)) + J— (2 sinh/^e^^ - l) e-^^* (/o(2a) - /2(22t)) (31 



Let us now investigate the asymptotic behaviour of the density for this (correlated) initial state (^9|). This is 
done as in the previous subsection. We begin with the massive, non-universal regime, where ht > \c\t ^ 1, 
and {b > c implies \z\ < 1) with u = LF'/\c\t ^ 1. In this regime, the density decays as p{t) — p(oo) = 

_il±IimE)f^:!:!l [e-'^sinh/3(l -fe-'^sinh/?) +0((|c|i)-i)]. 

4-\/ 7r|c|i 

On the other hand, in the critical regime (when e' = 0, 6 = c = e > 0, ct 3> 1 and z — 1, with u — 2L^/|c|i ^ 1), we 

. 1 — sinh /3(l+e~^ sinh /3 ) ^/-,// \_^/9\ 

have a power-law decay: p(t) — + C((c<) ) 

2V irct 

As before, in the case where there would be no pair-annihilation (i.e. h = — c = e' > 0), we have a power-law decay 
of the density: p{t) ^--^ 1 — -^^=. However in this case the asymptotic behavior is independent of p, and so is universal. 

Notice that in the critical regime, the density decays as a power-law : p{i) ^ as for uncorrelated initial cases 

(^,^^ . We also remark that the amplitude of the long-time behaviour of the density is non-universal and depends 
on the initial state through the parameter p (with = In -^)- 



We conclude that for the correlated states under consideration (|26||29|), initial correlations affect the dynamics 
through the amplitude of the density, which in absence of pair-creation is no longer universal. We have observed that 
when there are pairs created, without pair-annihilation, the density decays, for the two classes of correlated states 
( ^ , ^9|) , as an wmwersaZ power-law. 



D. Comparison with traditional free-fermion methods 



In this subsection we discuss and compare the approach devised in this work for the study of the DPAC model with 
ih.e free-fermion methods |J|,|,|l0|,|l]jl5|l|fl| . 

In most traditional free-fermionic methods, the stochastic Hamiltonian (|l|) is recasted, via a Jordan- Wigner trans- 
formation, in (free-)fermion form. The resulting quadratic (but generally neither hermitian nor normal) Hamiltonian 
is then reformulated in the Fourier space. As the resulting Hamiltonian is quadratic, the time dependence of the oper- 
ators is simple. Such an approach is well adapted for translationally invariant systems, and has been extensively used 
for the study of the DPA model (where e' = 0) with homogeneous initial density. Explicit results have in particular 
been obtained for initially empty and full lattice as well as for initial density p{Q) — 1/2 (see e.g. p|j2 pl] , p^ ). 

Another approach consists in diagonalizing the (free-)fermion version of the DPAC stochastic Hamiltonian through 
a Bogoliubov-type rotation, dealing, in Fourier space, with a diagonal quadratic form of the so called pseudo-fermions 
This method is particularly efficient for the DPAC model in the presence of pair-creation (e' > 0), but concrete 
results only have been obtained for lattice initially empty or completely filled 

Let us also mention that some exact results have also been obtained for particular non-uniform initial distribution 
for the DPA and DPAC models 

Here, to illustrate the approach devised in this work, in the translationally invariant situation, we have recovered 
known results for the density for both DPA and DPAC models (|l^,|l|,^ Hil!]) and have extended these results to the 
case of arbitrary initial density (|l^j2^,|2^,^,^^ . These results should, in principle, be also accessible by the traditional 
methods described above, but in the case of arbitrary initial density < p(0) < 1, the computations required are 
tedious. On the contrary, with the approach devised here, previous results are recovered in a simple and systematic 
manner. Although the long-time behaviour of the DPA does not depend on p(0), this not the case for the general 
DPAC in the presence of pair-creation, where the dynamics is no longer universal but depends on p(0) (through the 
amplitude and the power-law, see (|l^,^,^,|3[)). Therefore, when e > 0, e' > 0, the study of the long-time behaviour 
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of the density with respect of the initial density is relevant. An important advantage of the present method is the fact 
that it is based on an explicit generating function which allow the systematic computation of multi-point correlation 
functions. In addition, the method employed in this work provides explicit results which are beyond traditional 
methods: the results ( p^|Jl5|Jl^ ) are specific to this approach and are valid for arbitrary initial distribution (uniform, 
non-uniform, correlated, ...). Furthermore, this approach does not make use of the Jordan- Wigner transformation, 
and so we are not limited to deal with nearest- neighbour correlations (see (|29| - ^T1) ). Another advantage of the present 
approach is that it can be extended to the case of site-dependent DPAC model and will also be useful in the study of 
disordered version of this model 1^. In this sense the present approach is complementary to previous ones. 



VI. DENSITY AND NON-INSTANTANEOUS CORRELATION FUNCTIONS FOR DPACI MODEL 

In this section, we consider the model described by the stochastic Hamiltonian ijDPACi _ jjDPAC _|_ j^mput 
(according to (|2|)) where TJ^^^C? j^mput j-^^^yg been defined previously (in (|l|) and (^), respectively). From a 
theoretical point of view this model is interesting because it is one case where (non-trivial) nonequilibrium (the 
detailed balance is violated) steady-states can be computed exactly. 

To proceed with this study, following the original approach of Racz [|j , we first establish the equations of motion 
of the two-point correlation function of the dual counterpart of these models, i.e. TJ^^-^c (defined in (^)) and iJ™P"* 
(defined in (^), respectively. 

We obtain the equations of motion of the correlation functions as ^{(^m^m+i) = ^{['^m^m+i^ ^ iJ«"P«*]) 

For the translational invariant situation, we adopt the following notation (/ > 0): ri{t) = {(^m^^n+i) i^) = i^o'^Di^)- 
For the non-translational invariant situation, we define rm+i,m{t) = {^m+i^m){'^)- 

For the translational invariant system, the equation of motion of the two-point correlation functions is (c 7^ 0): 

= -2(& + 7;e(L-0)n(t) + c(n+i(0 + n_i(0), (32) 

with the boundary condition ro = 1, where Q[L — I) = 1 for I < L and Q{L — Z) = 0, for I = L. 

For the non-translationally invariant case, we consider the unbiased version of the DPAC model ((§) with h = h' >0, 
i.e. Lushnikov's model), and the equation of motion of the correlation functions then reads (I > 0): 

^rm+l,m = -2 [b + ^{2me{L - m) + I)] ri{t) + c{r,n+l+l,m{i) + 7'm+J-l,m(i) + r^n+Um + !(*)+ r^n+Lm-lit)) , (33) 

with the boundary condition ro,o = 1- 

i) We start with the translationally invariant problem (p2[). 

The steady state n corresponding to this situation obeys the difference-equation: 

2(6 + 7Ze(i - l))n - Z{fi+i + n_i) , (34) 

To solve (^), we make the Ansatz n = AJi^a {c/"f) ( where J^iuj) is the Bessel function of first kind) and take into 
account the following property of Bessel functions: Jv+i[u}) + Jv-i{u)) = ^J^{uj). 

Therefore: a = ^ — LSi^l. Taking into account the boundary condition ro = f(0) = 1, the constant A follows as: 

A ~ (^Jb{c/j)^ . The steady-state of the density for translationally invariant systems described by ( |3^ ) reads: r/ — 

— -^y^- . To solve the dynamical equation (^), we seek a solution in the form [|6|Jl^: ri{t) = n+X];/ 

With this Ansatz, ( p2|) reduces to the following difference equation: qr-i+i+'Zr-i-i = ^ [6 — cA/' — ^IQ{L — I)] qi'^i, 

which is solved as above, with an Ansatz: qii-i — BvJv-i+p^c/^). We thus find: (3 — ^ + 21 — V — L6i^l, where 
the spectrum A;/ is determined by the boundary condition ro = 1, which implies qii = 0, and imposes the condition: 

V^V)/7(^/^) = (35) 

The constants Bii are fixed by the initial condition: qi'-i = ri{0) ^ fi — BiiJ^^_~^ y~f+i-LSi l ^/^^ 

Finally, from the the spectrum of {Ac} defined by ( ^ ) and the Bii defined above , the solution of the equation (|3 
follows 
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(fc-cA,,)/7+;-L5,,i, 



(36) 



where the inverse of the relaxation time, say A* , is determined from the smallest zero of the Bessel function ( |35| ) and 
the steady-state has been obtained above. 

With the definition of the duahty transformation (^), as in (p^), we obtain for the density of particles p(t): 



l-ri(t) 



(37) 



where the steady-state is p{oo 



2 2 ' 



■^(i+V7)(g/7) 
~^Vt7(^7tO~ 



In the absence of pair-creation (e' = ^ 6 = c 



e > 



0), we recover the results of [|j (with b = c = e = 1). 

Following reference [||, we can now study the relaxation-time in the limit of a weak source (i.e. when 7^1) and 
analyze its interplay with the pair-creation term (described by (6 — |c|)/2). For the DPAC model, where e' > and 
7 = 0, we have seen in section V that the density (for translationally invariant systems) decays exponentially fast in 
time as ~ exp(— 2(6 — \c\)t). In ||] it has been shown that in absence of pair-creation (i.e. 6 = c = e > 0) the density 

of the DPAI model decays as ^ e^'^ Here we will study the case where both pair-creation (6 > c > 0) and the 
source term (7 > 0) are present. 

As we focus on the smallest value A* of the spectrum {|A;'|}, we consider 7^1. In this limit, we obtain an 
explicit relation which the {Ac} have to fulfill in terms of the zero's aj' of Airy function {Ai{ai') — 0), namely : 

(i-fiX y/^ ~ j ^'^^'-^ ^ = f ('^ith 1^1 > 1) and the a^/'s are real and negative. 



M 

(1-mAj)1/3 — \^2b2 j l"jh "i-^i- — b 

In the presence of pair-creation (i.e. when b > \c\), one has to consider the above equation , which should be solved 
for the (infinite) set of a^ 's. In so doing, we would obtain the spectrum {Aj}, which in turn provide the inverse of the 
relaxation-time A* = mm({|Aj|}). Here, we prefer to focus on 7^1, and /i ^ 1, i.e., we assume that the pair-creation 

1/3 



term is also "small" (e' <^ 1). Therefore, the spectrum obeys the following equation: A 



In the absence of pair-creation (i.e. /i = 1), we recover the result of |g] 
The inverse of relaxation time A*, follows (7 <C 1, e' <C 1 and fJ. ^ I) 



^3 - [2P 



3fi 
1/3 , 

\aj 



1 



A* = 



3c 



2b 



7 

262 



1/3 



(38) 



Thus the long-time behaviour {bt > \c\t ^ 1; 7 <C le' ^ 1 and /i < 1) of the density reads (see (pSh) : p{t) — p{oo) 



-e where B* corresponds to the term Bi' for which A;' = A*. 



ii) Let us now pass to the non-translationally invariant situation. In this case, we have to solve the linear differential- 
difference equation (^3|). To do this, we follow the same prodedure as above. The steady-state rm+;,m is inferred 
from the Ansatz fm+i.m. = A(Jm+i+a (c/7))^ and obtain: a — j- — ^ — LSm,L- Taking into account the boundary 



condition 



1 we find {I > 0): f„i+i. 



denote Xm'.i' simply by A. 

To solve the dynamical equation (|2 

y , 



_b_ 

27 



Hereafter for notational simplicity, we will 



-2Act 



we seek a solution of the form: rm+i,m{t) 



m'U 



'^m+l,m (^) 



and obtain the solution of 

-2Act 



m',l' 



,i'e 



+m+l-LS„ 



(c/7), 



(39) 



where the boundary condition qm'+v 



= determines the spectrum {A} through the condition: 



(40) 



The constants Am' v follow from the initial conditions. 



With (14, ^9|), we obtain the density of particles {m < L): 
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(41) 



m' ,1' 



with 



("'m+l)(oo) = 



1 - 



(42) 



iii) From the solution ( |4l] ) for the density in the non-translationaUy invariant situation, we can obtain the non- 
instantaneous two-point correlation functions of the (unbiased) DPACI model. 

With p-l<m<L, (n„+i(i)np(0)) - (n„+i(cx))) = - J2m' .1' %^e-2^=*Ji^ (c/7) ^^^(c/7) 
The only change from ( [4lD is in the B„i',i' which are fixed from the initial conditions. 

The spectrum {A} again follows from equation ( ^ ) and provides the inverse of relaxation time A*, as A* = 
min ({A„j}). 



With n 

b+2jn 



27 



LSm.Li in the limit j\n\ <ti 1, 

3/2 



e' ^ 1, with /I sa 1, the spectrum satisfy the equation: A,; 



c(2/i-l)3/2 



-1/31 



2/i?i 



^27-^ 1/3 



where the Uj < are again the zeros of the Airy function. 

It has to be emphasized that steady-states obtained for this DPACI model are nonequilibrium steady-states, in the 
sense that they violate the detailed-balance condition. In the DPAC model (where 7 = 0), the steady-states were 
expressed in terms of the quantity z ( with p(oo) = ^-j^) fi^'st introduced by Glauber [|l6[ to describe the steady-state 
of the time- dependent Ising model (with the prescription of recovering, as steady-state, the thermal equilibrium). Thus, 
in this sense the steady-state of the DPAC model is reminiscent of the states in equilibrium statistical mechanics. On 
the contrary, the steady-states of the DPACI model are genuine nonequilibrium steady- states. 



VII. CONNECTION WITH THE COAGULATION MODEL 



In this section we translate the results obtained so far for the DPAC-model to another, free-fermionic, model 
of stochastic hard-core particles on a lattice, namely the diffusion-limited coagulation (DC) model. In this model 
particles can jump from a site r to an adjacent site (provided it was vacant) with rates h {% + A — > A + 0) and h' 
{A-\-% — > -I- ^). Two adjacent particles at site r and r — 1 can coagulate with rate h {A + A — > A + %). Two 
adjacent particles at site r and r — 1 can also coagulate with rate h' {A-\- A — * % A): 

It is known that (for initially uncorrelated systems) there is an exact duality between DPA-model (where there no 
pair-creation: e' = 0) and the coagulation model (see e.g. P, pl] , p^p^ ] and references therein). Here we briefly recall 
the transformation connecting these two stochastic models and translate the results of the DPA-model. 

Applying the simikrii?^transformation ||,lll|,|l|,p2| H = BHB^^ ; B = B ® B ® . . . ® B = B'^^ , where 



B — 2 j '' (stoc/iasti!c-)IIamiltonian H of the DPA-model becomes the (stoc/iasfi!c-)IIamiltonian H of the 

coagulation- model. 

For the coagulation model (with initial state |P'(0))), we can express the correlation functions with help of cor- 
relation function of the DPA model. The density of the DC then reads: (rij(t))|'^9'('Q')^ — 2{nj{t))^f/^p,^g,^y For the 

two-point correlations functions, we have:(nj+.r(i)nj(t))|'^9'(o)) = 4(nj+r(^)7^j(i))e^L"|p/(o)) ^'^'^ ('^j+r(0"j(0))jp'?(o)) — 
2{nj+r{t)nj (0))g^'i'|p,(Q-,^ 

As an illustration, let us no w p rovide, the density of the DC model from the expressions obtained for the DPA- 
model (where e' = 0) as (see ([14| )). The density of the translationally invariant DC model with uncorrelated initial 
state |P'(0)) = (1 — p{0) p{0)y (of initial density p), is related to the unbiased version of the DPA model with initial 
density Thus, in this situation, with ( [T^ ) , we obtain for the density of DC model: 

p^C(^) ^ ^-2(h+h')t I ^^(2(/i -I- h')t) + h{2{h + h')t) - ^(1 - p{Q)Y (/,_i(2(/i -I- h')t) - Is+i{2{h + h')t)) \ (43) 

I s>Q J 

As for the DPA-model, for translationally invariant systems, the instantaneous correlation functions do not depend 
on the bias h — h' . We could also extend in a similar manner all the results obtained in the previous sections V to the 
DC model. The discussion of the long-time behaviour of p^'~''{t) follows directly from that of the DPA model given 
in subsection V.B. 
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VIII. SUMMARY AND CONCLUSION 



This work provides a complete and unified solution of some related one-dimensional nonequilibrium models. 

In a first part we have solved the diffusion-limited pair- annihilation- creation (DPAC) model (under the free-fermionic 
constraint) by a novel approach, alternative and complementary to the traditional free-fermion methods. Combining 
domain- wall duality with a recently developped technique |^ , we were able to compute the generating function 
of the dual of the DPAC model (section IV), i.e. to obtain correlation functions (section V). In this paper, we take 
advantage of this approach to recover results previously obtained (see (|l9| , pl| , p2|j2q )) by traditional methods, extend 
these (see(|l^,|^|3 24 25|j27|)) and derive original ones (see(|4|,^ 16 30|j31|)), in a simple manner. The method is very 
promising for site-dependent and disordered systems j^. We focused on the density, the staggered- current, and the two- 
point correlation functions (instantaneous and non- instantaneous) (see ( |l4| , ^ |l6 ) ) . In particular we studied explicitly 



the density (in absence as well as in presence of pair-creation) for translationally invariant systems with arbitrary initial 
density and have discussed the asymptotic behaviour of these quantities in different regimes. We obtained explicit 
expression of the staggered- current in arbitrary translationally- invariant (uncorrelated but random) initial conditions 
(see (^Hj)). We investigated the densities for two classes of correlated (but tranlationally invariant) systems. We have 
computed explicitly the non-instantaneous two-point correlation functions for initially translationally invariant and 
uncorrelated (but random) states: these quantities depend on an eventual bias (contrary to instantaneous correlation 
functions). As a by-product we provided the solution of a biased generalization of Glauber's model, which is the dual 
counterpart of the DPAC model. 

In a second part we studied the density and non-instantaneous two-point correlation functions of the diffusion- 
limited pair-annihilation-creation model in the presence of source of particles (DPACI model) . Via the domain- wall 
duality transformation, we solved exactly the equations of motion of the correlation functions of the dual model. 
We obtained the exact density and steady-states for translationally invariant as well as non translationally-invariant 
systems. We also computed the non-instantaneous two-point correlation functions. In particular for the translationally 
invariant situation in the absence of the pair-creation term, we recovered, the density first computed by Racz |^ . We 
extended these results to the case where both pair-creation term (e' > 0) and source term were present (7 > 0) and 
obtained the relaxation-time for the case when both 7 <C 1 and e' <C 1 (see (38|)). 

In the last part of this work, we have shown how to translate the resul 
diffusion-limited coagulation (DC) model. 



ts obtained for the DPA model to the 
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